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We study the numerical propagation of waves through future null infinity in a conformally com- 
pactified spacetime. We introduce an artificial cosmological constant, which allows us some control 
over the causal structure near null infinity. We exploit this freedom to ensure that all light cones 
are tilted outward in a region near null infinity, which allows us to impose excision- style boundary 
conditions in our finite difference code. In this preliminary study we consider electromagnetic waves 
propagating in a static, conformally compactified spacetime. 



I. INTRODUCTION 

The motivation for the present work is the matching of 
a computational model of a source of gravitational waves 
(such as a compact object binary) to the cosmological 
environment through which any emitted signal travels to 
a detector. The transmission phase is well understood 
and is essentially the propagation of a linearized grav- 
itational wave through a background cosmology which 
includes such effects as a cosmological redshift and grav- 
itational lensing due to intervening masses, plus the ki- 
netic effects of the detector's motions. For these purposes 
the geometrical optics approximation is often adequate. 
Although implementing this transmission model in ob- 
servatory data analysis demands skill, it is simpler than 
the electromagnetic case where interactions with dust, 
gas and plasmas en route must also be considered, fn 
this large scale background the source can be taken to 
be a point source of linearized gravitational waves hav- 
ing some specific antenna patterns and wave forms. The 
more unsettled question then is how to extract from a 
detailed model of a wave generation system the appro- 
priate linearized wave. The most common assumption is 
that the wave generator can be modeled as existing in an 
asymptotically flat spacetime where waves propagating 
out to future null infinity (scri-f or ^+) will be described 
there by linearized gravitation, originating from a point 
source in the cosmological context. This is clearly jus- 
tified by the huge difference in scales between the size 
of the wave generator and the radius of curvature of its 
surrounding spacetime (the cosmos or a galaxy or stellar 
cluster) . The question then is how to efficiently calculate 
and extract the waveforms at 

In the usual Cauchy evolutions on asymptotically flat 
spacelike hypersurfaces, the waves must be recognized at 
large distances from the source, possibly a long time after 



the source motions have ceased violent activity. One way 
to overcome this time lag is through the use of a retarded 
time coordinate. This works well in spherical symmetry 
Q> Q J but is poorly defined otherwise. The Pittsburgh 
group has used such characteristic evolution at larger dis- 
tances while matching to more common Cauchy evolution 
in the interior The use of hyperboloidal time slices 
(some generalization of hyperboloids in Minkowski space- 
time) appears attractive as a smooth transition from con- 
ventional time slices in central regions to asympotically 
null spacelike hypersurfaces toward infinity. Were such 
a slicing to be combined with a straightforward Cauchy 
evolution, the Courant time step condition would defeat 
the advantages of approximating retarded time, as the 
nearly infinite coordinate speed of light toward infinity 
would require infinitesimal time steps (since the Courant 
condition generally requires that Ai/Ao; < 1/c for nu- 
merical stability). Thus hyperboloidal slicings effectively 
require a conformal compactification which makes the co- 
ordinate speed of light finite again 

In addition to the use of hyperboloidal time slices and 
conformal compactification, we suggest two additional 
tools to ameliorate the extraction of outgoing waves in 
numerical evolutions. One is the use of artificial cosmol- 
ogy to convert the null hypersurface into a spacelike 
hypersurface in the future of the initial Cauchy hypersur- 
face as proposed in . This procedure fattens J^'^ into a 
shell of spacelike hypersurfaces by changing the value of 
the cosmological constant, plausible for the present uni- 
verse, by many orders of magnitude so that the de Sit- 
ter horizon is located, not at cosmological distances, but 
merely well beyond the region where the wave generating 
mechanism can be influenced by it. In this sense it is sim- 
ilar to the computational tool of artificial viscosity which 
fattens a shock from its physical thickness on the order of 
the molecular mean free path into a much larger length 
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which is still small in the range of scales that otherwise 
are met in the physical problem at hand. The second 
additional tool is the tilting of the null cones governing 
wave propagation in ways that ease the computational 
behavior in the unphysical ('das AH') region beyond the 
de Sitter horizon and beyond . 

The German noun "das All" refers to the Universe, 
but sometimes in a sense unfamiliar in English where 
it translates better as "outer space." A spacecraft, af- 
ter launch, might be described as entering das All. In 
the computations reported here, it is precisely this outer 
space — the physical details far from a source of gravita- 
tional (or in our case electromagnetic) waves — which is 
excised from the computational model in analogy to the 
way the internal mysteries/singularities of a black hole 
are often excised in computational models. We develop 
and implement the proposals made earlier in Refs. 0, Q 
and suggested above to obtain wave forms promptly at 
infinity by using hyperboloidal time slices in a Cauchy 
evolution, to avoid infinitesimal time steps through con- 
formal compactification, to convert outer boundary con- 
ditions to excision boundary conditions by using an ar- 
tificial cosmological constant, and to defend the physical 
domain against computational boundary errors by tilt- 
ing the light cones further outward in the region beyond 
the outer horizon. The aim of this paper is to explore 
the application of this set of boundary treatment tools 
in a low cost application in order to see whether there 
are unexpected impediments to its application, and to 
get guidance for the best approaches for its application 
to larger problems. 

The use of hyperboloidal slicings in numerical relativ- 
ity has been carefully studied (e.g. Refs. iS, iTi, |8|| ) , and 
there are existing numerical results generated on com- 
pactified domains (e.g. Ref. 'o']). This paper contributes 
to this literature by combining these two techniques in 
a spacetime with an artificial cosmological structure de- 
signed to make it easier to construct boundary condi- 
tions at J^+j and by studying the effectiveness of mod- 
ifying the light cone structure in the region of the ana- 
lytically continued spacetime exterior to J^"*" to prevent 
non-physical, incoming radiation from piling up on the 
non-physical side of J''^ . That said, our approach to 
the problem does not require special modifications to the 
form of the differential equations solved to model the 
physical problem. It should be compatible with existing 
numerical codes desgined to solve the initial problem, de- 
scribed by some variant of the the Einstein equations in 
an Arnowitt-Deser-Misner 3-1-1 decomposition 1^. 

Our work is also complementary to community wide 
efforts to use mesh refinement techniques in numerical 
relativity. Several groups have recently and successfully 
incorporated mesh refinement techni ques into their Ein- 
stein solvers (i.e. Refs. H [H IH El 111 B E E El 
0). This technology allows groups to push their com- 
putational boundaries to very large coordinate distances, 
reducing the effects of incoming radiation and constraint 
violating modes introduced by applying boundary condi- 



tions at a finite distance. 

While mesh refinement has greatly increased the power 
of modern codes, it also has limitations. It is now fairly 
easy to surround a computational domain with coarser 
and coarser grids, but at some point the coarse grids 
cannot resolve a wavelength, see, e.g. Refs. E3 a-^d |2l 
Section 4.3] . They can, however, insulate the active, wave 
generation regions of the computational domain from 
boundary problems for a limited period of time. During 
this period waveforms have been successfully extracted 
from inner regions as small as 2.4 wavelengths from the 
center E3- The hyperboloidal slicings considered here 
and elsewhere have the advantage that they asymptote 
to the outgoing light cones in the spacetime. Outgoing 
waves, therefore, appear asymptotically constant and re- 
quire very few points to resolve them all the way to J"^ . 

The rest of the paper is organized as follows: In Sec.lTTI 
we briefly summarize the formalism we use in our code, 
both for the background spacetime metric and for the 
Maxwell equations themselves. In Sec. IIIII we describe 
the numerical methods used to solve the equations. We 
summarize the numerical results and discuss the impli- 
cations for future work in Sec. IIVI Additional detail is 
provided in the Appendix. 



II. FORMALISM 

A. Spacetime Metric 

As in Ref. the background metric for this study will 
be the de Sitter spacetime 

ds^ = ~dT'^+dX'^+dY^ + dZ^ + {R^/L^){dT^dRf (1) 

with an artificial cosmological constant A = -|-3/L^. We 
are interested in the limit when this becomes 

Minkowski spacetime, and the cosmological constant is 
used to make boundary conditions numerically simpler 
(we hope) at the de Sitter horizon than they would be 
at fiat spacetime's The de Sitter horizon is the null 
hypersurface R = L in this metric. We introduce the 
coordinate changes 
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s 1 - r2/4 

as used in Refs. where R^ = X^+Y^ + Z"^ = X'X\ 

= a;^+T/^ + z^ = cc'a;^ and s is a constant scaling factor 
whose geometrical significance can be read from Eq. (O 
below. This brings in to r = 2 in the Minkowski case 
and in the de Sitter case makes this a spacelike hyper- 
surface beyond the de Sitter horizon which is part of de 
Sitter spacetime's future causal boundary which is also 
called . The hypersurfaces of constant u are then 
hyperboloids 



[r-s(u~l)]2-i?2 = s2 



(3) 
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in Minkowski spacetime and are also asymptotically null 
spacelike hypersurfaces in the de Sitter modification. 
These hyperboloids in Minkowski spacetime have con- 
stant curvature of radius s which is therefore the scale of 
the region in which they are approximately flat. Also the 
time interval AT at the spatial origin between the u slice 
on which a light signal is emitted there and the u slice 
on which it reaches is Au = 1 or AT = s. Thus s 
may also be called the scri-delay. This coordinate change 
leads to a metric which is singular at r — 2, but only in 
a conformal factor s^/g^ = — which does 

not appear in the Maxwell equations. Thus we can drop 
this conformal factor and our test problem is to solve the 
Maxwell equations in the resultant metric which is of the 
form 



-a^dt^ + 7y {dx' + P'dt) {dx' + P^dt) (4) 



with t — u. For use in later modifications we include a 
function W{r) in the definition of the metric functions in 
0). For now we choose W = and then the following 
equations give the analytically continued and conformally 
regulated de Sitter metric described above: 
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To see the behavior of the light cones for this regulated 
metric, we calculate the coordinate speed of light dr/du 
in radial null directions. This gives 



~ 1 + (sVL2)52 

For the de Sitter values above this reduces to 



and 



Vout = (1 + \rf 



l + (s2/i^)'S'2 



(7) 



(8) 



(9) 



Note that the constant r hypersurfaces are never space- 
like when s/L = since then the inward speed of light 
is Vin < 0. When the de Sitter option is used one finds 
Vin > 0, i.e., both the inward and the outward sides of 
the lightcone point toward increasing r, for an interval 
around r = 2. 



In evolution codes which permit setting the computa- 
tional boundary on a coordinate sphere (such as those 
described in [2ll IH, 111]) the outer boundary could 
be set at r = 2 or slightly beyond in the metric described 
above (or ones asymptotically similar). However, if a cu- 
bic boundary enclosing ^+ is to be used, there are addi- 
tional behaviors beyond that need attention. For ex- 
ample, if the artificial de Sitter horizon is set at L = 10s, 
then the r = constant hypersurfaces are spacelike only 
in the narrow interval of (approximately) r = 1.81 to 
r — 2.21. But a boundary cube enclosing r = 2 has cor- 
ners at r = 2^/3 « 3.46. Thus on most of such a cubic 
boundary the light cones (from Eq. ^) would allow in- 
wardly propagating waves, requiring careful attention to 
the boundary conditions which one hopes to avoid. 

Fig.Qlbelow shows the coordinate speeds for light rays 
on the inward side of the light cone Uin, on the outward 
side of the light cone Wout, and for a timelike center of 
the light cone Wccntcr = ~/3 normal to the time slices 
of constant u. This allows light rays to move inward 




FIG. 1: Coordinate speeds for the de Sitter lightcone: inward 
(uin, blue, lowest curve), normal to the spacelike time slices 
(— /3, green, middle curve), and outward (uout, red, highest 
curve). For the Minkowski metric the graph looks very simi- 
lar, except that Vin is there never positive. 

(win < 0) in some of the region beyond J^'^ within a 
cubical computational grid. 

In this case we suggest that the metric be modified 
outside r = 2 so that all r = constant hypersurfaces be- 
yond r = 2 are spacelike — i.e., so that > in the 
entire computational domain beyond r — 2. Analytically 
this would have no effect on the physical solution inside 
r = 2 provided all fields propagate causally. This mod- 
ification is also straightforward to do in the test case of 
a Maxwell field with a fixed background metric. If these 
methods were to be developed for the Einstein equations 
(without using a spherical boundary such as r = 2) it 
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would necessitate distinguishing, in the Einstein equa- 
tions beyond r = 2, the g^i^ metric being evolved from 
the appearances of the metric as coefficients of the prin- 
ciple derivatives which determine the causal relations in 
the evolution algorithm. 

In the Minkowski case s/L = (no artificial cosmol- 
ogy) it is not possible to make a change in the metric 
only beyond r = 2 which avoids inward pointing null rays 
just beyond r — 2 — in that case one has, from the r < 2 
formula above for that {(f/dr'^)vi^ = — ^ at r = 2 so 
that Via must have a local maximum value of zero at that 
point and is necessarily negative (inward) in a neighbor- 
hood of that point, including some r > 2 points. For 
the de Sitter case there is again a maximum of Win near 
r = 2, but its value there is positive so it can be smoothly 
turned upward (to avoid zero or negative values) in suit- 
able intervals beyond the de Sitter horizon (and beyond 
r — 2 for small values of s/L). 

In Eqs. © and assuming s/L > 0, we choose to 
take 

(r - 2)"^ 

W{r)^e{r-2) ^ '- 2 (10) 

to be a smooth step function that modifies the metric 
beyond r = 2 in a way that keeps the "incoming" coor- 
dinate speed of light Vin = dr/du positive (that is outgo- 
ing) when r > 2. The function is a C"^ modification 
of Q[r — 2), the Heaviside step function that jumps from 
to 1 at r = 2. Fig. |21 below shows the modified light 
cone directions in this case. 

Note that for r < 2, where W — Q, this metric has 
the surprising property that the outgoing speed of light 
Vout = (1 + ''/2)^ is independent of the cosmological con- 
stant parameter L. Although in general the presence of 
a cosmological constant alters the causal structure of a 
spacetime, and indeed it does modify the ingoing radial 
speed of light in the present example, it does not change 
the outgoing radial speed of light from that of the origi- 
nal Minkowski spacetime. This fact will prove useful for 
numerical diagnostics. 

This metric surgery has an additional side effect worth 
noting. For this modified de Sitter metric using Eq. l(Tn|l 
one finds very low inward light velocities in the region 
just outside as seen in Fig. 13 This can lead to fields 
propagating along the inward side of the lightcone spend- 
ing very long periods of time in the region r « 2.3 as 
illustrated by the world line of a light ray in Fig. 0] 

B. Maxwell Equations 

We take the Maxwell equations in the form given in 
the DeserFest paper |^. The quantities involved are all 
various components of the usual 4-dimensional Maxwell 
fields F^^ and Tf"" = y^gt'^'g^'^Fap. In particular, for 
any 3-1-1 metric of the form Eq. Q we define, 

jro^ = V' , T'^ = [ijk]Hk (11a) 




FIG. 2: Coordinate speeds for the modified de Sitter (^ean.llOfl 
lightcone: inward light rays, {vh-i, blue, lowest curve), motion 
normal to the spacelike time slices (— /9, green, middle curve), 
and outward light rays («out, red, highest curve). 

v_in (r) 




FIG. 3: Coordinate speeds inward (t^in) for fight rays in the 
modified (egn. ITUll de Sitter metric. 

and 

F.o^E, , Fj^[ijk]B'' . (lib) 

Here [ijk] is the completely antisymmetric symbol [ijk\ — 
0,±1 with [123] = +1. The fundamental fields in the 
formulation are and I?*. The constraints read 

d.,B' = d{D' . (12) 

We need only take regular partial (rather than covariant) 
derivatives in our formulation of Maxwell's equations be- 
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dr^^^ilin (r) 




50- 



O j i' 

1.8 2 2.2 2.4 2.6 2.8 

r (u) 

FIG. 4: World line of a light ray on the inward side of the 
light cone in the modified de Sitter metric. 

cause of judicious use of tensor densities rather than ten- 
sors (cf. Exercises 22.8 and 22.9 in [l^). The resulting 
fields in Eqs. (|12|) are each three dimensional vector den- 
sities so this equation is three dimensionally covariant. 
In terms of these primary fields, two auxiliary fields Ei 
and Hi are computed by the formulas 

E^ = {a/^UjV^ + [ijk]f3^B'' (13a) 

and 

H, = {a/^)j,^B' -[ijk]P^V'' . (13b) 
The evolution equations are then 

doB' = -[ijk]djEk (14a) 

and 

doV' = [ijk]djHk . (14b) 

Note that the metric does not appear in Eq. (|12|l or 
Eqs. H14(l . and that the metric quantities which do ap- 
pear in Eqs. H13|) are conformally invariant. This is a first 
order partial differential equation system in flux conser- 
vation form. 



C. Initial Data 

For initial data and some testing purposes we have 
used the solution of the fiat spacetime Maxwell equations 
employed by Knapp et al. |2g and subsequently by Fiske 
[27| . This solution is a wave pulse, which propagates 
smoothly from through the origin and out to 



It consists of a toroidal T> field and a poloidal B field 
generated from a vector potential, described in Ref. 5] , 

A ^ A,dX' ^ f sin^ ed(t> (15) 

where 

sin^ ed(j) = -X: (XdY = YdX) = 4" (xdy - ydx) (16) 
and 

/ = - 2AL/) er^^' - + 2XV^ e-^^' (17) 

where the constant A parameterizes the pulse width. In 
these equations U = T ~ R and V = T + R are the 
Minkowski retarded and advanced times, which are then 
expressed in terms of the compactified coordinates u and 
r as U/s = u-2 + 2/{l + \r) and V/s = u-2 + 2/ {l-\r). 
Finally, B"^ and T)^ are calculated by differentiating this 
A field. The term containing e~^^ can be smoothly 
omitted when r > 2 since all its derivatives vanish at 
r = 2, but the term with e~^^ must be retained to 
preserve the constraints, and we thus have nonzero initial 
data in the unphysical region r > 2, i.e. beyond J^~^ ^ 
where initial data are not physically meaningful. As the 
initial pulse is centered at r = 0, we would ideally like to 
make it sufficiently narrow so that its amplitude would 
rapidly taper off and become negligible beyond r > 2. 
However, we did not achieve enough resolution to reduce 
the pulse width to make these unphysical data beyond 
r > 2 insignificant. 

For all of the simulations described in this paper, we 
choose A = 1 (or, equivalently, we measured distances in 
the uncompactified space in units of A^^/^). The pulse 
width is then A~^/^. We also always choose s = 1 to set 
the size of the region near the origin where the hyper- 
boloidal slices are somewhat fiat. Our results reported 
below used either s/L = for the Minkowski case, or 
s/L = 0.1 for the de Sitter and modified [Eq.HUl de Sit- 
ter cases. Thus the r = constant hypersurfaces were not 
spacelike for the de Sitter case when r > 2.21 

III. NUMERICS 

A key point in our approach to this problem is that we 
do not require special numerical techniques to handle the 
compactified spacetime at or within J^^. At all interior 
points, we use second-order, finite-difference approxima- 
tions to all derivatives. At the computational boundary 
beyond J^^, we use a second-order, one-sided stencil for 
normal derivatives. We integrate forward in time using 
the iterated Crank-Nicholson method commonly used in 
numerical relativity codes [28| . 

Our decision to use one-sided differencing at the com- 
putational boundary was (in our view) the simplest thing 
to do. It turns out, however, to be equivalent to the 
black hole excision boundary conditions first studied in 
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FIG. 5: The L2 norms (over the region r < 2) of the con- 
straints plotted vs. time. The higher resolution data is multi- 
plied by the factor appropriate to demonstrate second order 
convergence. The initial values arise from B and "D fields 
which are exact analytic solutions of the constraints, eval- 
uated at the grid points and then differenced to form the 
constraints. 



Ref. [23 ■ To be precise, along a normal to the boundary, 
the composition of third order extrapolation to a cell just 
outside the boundary (as used for excision) with second- 
order center differencing at a cell just inside the boundary 
gives the same stencil as second-order one-sided differenc- 
ing at the cell just inside the boundary. This is consistent 
with our mental picture of excising a region exterior to 
J^"*", and reflects the similarity between (at which 
all characteristics point outward) and an apparent hori- 
zon of a black hole (at which all characteristics point in- 
ward). As with black hole excision, it is critical for stable 
and accurate simulations that all light cones at the exci- 
sion boundary tilt towards the excised region, hence the 
need for the beyond-^^ metric modification described 
in Sec. Iff Al Without these Eq. HIO() modifications, the 
light cones would point outward only in a neighborhood 
of which did not extend out to the computational 
boundary of the cubic grid enclosing . 

We have verified that our code converges at second 
order accuracy throughout the physical domain. We paid 
special attention to the region in the vicinity of J^"*" , and 
we are satisfied that the code converges there as well. 
Fig. explicitly shows a two-point convergence test on 
the constraints over the physical region r < 2, which are 
seen to be second-order convergent to zero and nearly 
constant. 



IV. RESULTS 

When evolving with the original, conformal, de Sitter 
metric given in Eq. large errors eventually appear 
near the computational grid boundary, grow exponen- 
tially in time, and propagate inside J^+. This problem 
was slightly worse when the cosmological constant was 
omitted. These cataclysmic errors originated at the grid 




FIG. 6: A spacetime graph of r^B^ . Notice that the wave 
passes through J^"*", located at r = 2 in our coordinate system. 
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FIG. 7: The time evolution of B^ at a point on J^^ . We show 
both a numerical solution generated by our code using the 
modified de Sitter (engineered lightcones beyond J^^) and the 
analytic solution with zero cosmological constant. The close 
agreement of the lines suggests that we are able to accurately 
capture the behavior of the fields as they pass through . 
Other numerical evolutions based on the original (compact- 
ified, analytically continued) Minkowski metric (s/L = 0.0) 
and the similar de Sitter metric {s/L = 0.1), for which the 
grid boundaries were not spacelike, fail after t = u — 4. 



boundary which in these cases was not a spacelike hyper- 
surface for which excisions (at apparent horizons) were 
designed. As anticipated in Sec. Ill Al we proceeded to ar- 
tificially tilt the lightcones outwards beyond to make 
the grid boundaries spacelike. 

Fig. shows a spacetime view of the z component of 
the magnetic field for an evolution using the modified 
fean. llO|l de Sitter metric. In order to make clear that the 
wave is cleanly passing through ^+ through the course 
of the evolution, we have multiplied the field by , mak- 
ing the scale of the field comparable at the origin and 
at J^'^ . This can be seen in a more local way in Fig. |7| 
which shows the time evolution of at a single point 
on J^"*", for four cases. One is the analytic value of the 
field (in extended Minkowski spacetime) for comparison. 
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The two cases which failed to maintain vanishing fields 
in the wake of the outgoing pulse correspond to the two 
(Minkowski and de Sitter) metrics for which the modifi- 
cations needed to make the grid boundaries beyond 
spacelike [Eq. IjlOl) ] were not implemented. 

Note that the analytic solution is with zero cosmolog- 
ical constant, whereas the numerically evolved solution 
is with finite cosmological constant, and yet the wave- 
forms coincide. This phenomenon can be understood as 
follows. Because, as pointed out in Sec. Ill Al the outgo- 
ing speed of light is the same in either case, the phase of 
outgoing waves can be said to be independent of a cos- 
mological constant here. Further, since by construction 
the initial wave pulse is identical in our Minkowski and 
de Sitter simulations, by conservation of energy the wave 
amplitudes at are expected to be commensurate. 

These results show that many desired aims of this ap- 
proach are achieved. Waves propagate smoothly through 
^+ and no difhculties appear to arise in the neighbor- 
hood of J'^ The compactified hyperboloidal slices allow 
waves to appear at J^^ a very short computational time 
after they originate in the central region. No modifi- 
cations were needed to integration algorithms designed 
for conventional spacelike Cauchy slicings. The excision 
at a cubic computational boundary, however, only main- 
tained acceptable behavior for a time several times the 
pulse width, and ultimately led to behaviors in the un- 
physical region beyond ^+ which became intolerable. 
An example is given in Fig I A2I where the constraint vio- 
lations beyond ^+ are seen eventually to increase expo- 
nentially. In addition, by using analytically constructed 



initial values satisfying the constraints, we had to accept 
nonzero initial conditions in the unphysical region r > 2 
beyond J^+, and the subsequent evolution of these fields 
complicated the computation beyond J'^ . 

We believe that further developments of Cauchy evo- 
lutions with compactified hyperboloidal slicings and ar- 
tificial cosmology should be done with the computa- 
tional/excision boundary a spherical hypersurface at or 
modestly beyond . The means for doing this have 
been developed, e.g. [Ullllllllll. Then the possibility 
remains open that well posed wave equations could run 
indefinitely. As explained in the supplementary material 
(Appendix), here we did not use the differential equa- 
tions from Sec. VI of [s^l as our lightcone engineering 
(see Sec. Ill A|) beyond J'^ makes the constraints there 
increase exponentially even from analytic arguments. A 
spherical boundary could be spacelike at r > 2 without 
the need for such artificial tilting of the lightcones, and 
thus allow the use of the formulation Equations 26] 
where constraints evolve causally. 
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APPENDIX A: AUXILIARY MATERIAL 

This section is mainly a collection of graphs which il- 
lustrate properties that are merely asserted in the main 
article. 



1. Constraints 

Even when the modified de Sitter metric is used, the 
solution in the unphysical region deteriorates as shown 
in Figs. EU and US 

In spite of these problems, the field at ^+ remains zero 
after the wave pulse passes (Fig. [Tj) for a moderately long 
time, which can be lengthened by using higher resolution. 
See Fig. IA3I Better control of the constraints should be 
possible using the equation system 30, Equations 26] 
where the constraints propagate inside the lightcone and 
damp at a rate —di(3^. However this is damping only for 
the Minkowski or de Sitter metrics as seen in Fig. IA4I 
and becomes anti-damping for our modification of the de 
Sitter metric (see Fig. IA5I ) which makes our cubic grid 
boundaries spacelike. 

2. Initial conditions beyond J^"*". 

A second reason which makes the use of the full cubic 
grid (rather than only a region inside a spherical bound- 
ary such asr = 2orr = 2.1) inappropriate for numerical 
computation is the need for initial values beyond 



that is beyond r = 2, which is i? = oo. 

With a spherical excision boundary and causal prop- 
agation of all modes, there is less room for bad initial 
conditions to intrude, and better reason to think that 
even poor initial conditions would be quickly fiushed out 
of the computational domain. We have simply used as 
initial conditions the formulas for B and T) as computed 
from equations H15(l through (|17|l extended analytically 
by the coordinate transformation from TSY Z to uxyz 
of equations Q . As a solution of the Maxwell equations 
in the Minkowski metric, these formulae give substan- 
tial activity in the r > 2 region, as seen in Fig lA6l In 
that Figure one solution has been modified by smoothly 
deleting the terms involving exp —\V^ for r > 2 which is 
possible since V/ s = u — 2 + 2/(1 — ^r) becomes infinite 
at r = 2. However the pulse seen moving inward from 
beyond the grid will have depended upon initial condi- 
tions far outside r — 3.64 = 2.1^/3 which was the limit of 
our grid. Our initial data, however, are not evolved us- 
ing the Minkowski metric, but usually with the de Sitter 
or modified de Sitter metrics, where the decomposition 
into ingoing and outgoing parts of the initial conditions 
will differ seriously from the Minkowski case, especially 
well beyond Thus we can understand that differ- 

ent, but equally visible, activity could occur in the r > 2 
region which would be actually solving the given differen- 
tial equations but be unrelated to the physical activity in 
the region r < 2. This appears to be occuring in the nu- 
merical results (modified de Sitter for which no analytic 
solution is available) shown in Fig. IA7I 
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FIG. Al; Samples of the time evolution of the constraint violation Idivl?! at 2; = within our comptational cube \x'^\ < 2.1 
whose edges are at r = 2.97 . 



r>2 
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FIG. A2: The time evolution of an L2 norm of the constraint 
violations over the unphysical region r > 2 beyond J^"'" within 
our computational cube \x^\ < 2.1 whose corners are at r = 
3.64 . 
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FIG. A3: A continuation of the test wave pulse to later times 
at J*^+. 




FIG. A4: The damping factor for constraints in the de Sitter 
metric. 
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FIG. A5: The damping factor for constraints in the modi- 
fied de Sitter metric. Where div/3 is positive, the constraints 
should grow exponentially at the plotted e- folding rate. 
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FIG. A6: For the analytic solution in Minkowski spacetime: 
The magnetic field at z = after the pulse from Fig |7| has 
exited the computational grid. 
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FIG. A7: A field component plotted along a line z — 0,x — 1 parallel to the y-axis. At the grid limit y — 2.1 one has r = 2.33, 
while J^^ is met at y = 1.73. Inside the field values appear to be convergent to zero after it = 4, while in the region beyond 
y = 1.85 or r = 2.1 a slowly moving pulse appears to remain independent of the selected resolution. 



